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Abstract 

We study soliton solutions of the Nicole model - a non-linear four- 
dimensional field theory consisting of the CP^ Lagrangian density to 
the non-integer power | - using an ansatz within toroidal coordinates, 
which is indicated by the conformal symmetry of the static equations 
of motion. We calculate the soliton energies numerically and find that 
they grow linearly with the topological charge (Hopf index). Further 
we prove this behaviour to hold exactly for the ansatz. On the other 
hand, for the full three-dimensional system without symmetry reduc- 
tion we prove a sub-linear upper bound, analogously to the case of the 
Faddeev-Niemi model. It follows that symmetric solitons cannot be 
true minimizers of the energy for sufficiently large Hopf index, again 
in analogy to the Faddeev-Niemi model. 
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1 Introduction 



In the last few years there has been rising interest in non-hnear field the- 
ories which allow for the existence of ringlike, or, more generally, knotlike 
solitons. On the one hand, this interest is due to the fact that there may 
exist physical applications for such models, as is the case, for instance, for the 
Faddeev-Niemi model, which finds some applications both as a candidate for 
a low-energy effective theory for Yang- Mills theory and in condensed matter 
physics. On the other hand, the rising interest is related to the advance of 
more powerful computer facilities which allow for reliable numerical calcula- 
tions of those solitons in cases when an analytical solution is not available 
(which happens quite often). In addition, there is some intrinsical mathe- 
matical interest in theories with knot solitons. 

In the simplest case the field of the theory describes a map from one- 
point compactified three-dimensional space IRq to the two-sphere S'^. JRq 
is topologically equivalent to the three-sphere S^, therefore such maps are 
characterized by the third homotopy group of the target space S'^, which 
is nontrivial, 773(6'^) =2. As a consequence, fields which describe maps 
IRq — > fall into different homotopy classes, and a soliton is a field config- 
uration which minimizes a given energy functional within a fixed homotopy 
class. The topological index characterizing the homotopy class is called Hopf 
index, the corresponding map is called a Hopf map, and the minimizers are 
sometimes called Hopf solitons. For details on the Hopf map we refer to 
Appendix A. 

The probably best-known theory which allows for Hopf solitons is the 
Faddeev-Niemi model ([1], [2]) with Lagrangian density 

£fn — ^2 — (1) 
where A is a dimensionful coupling constant, £2 is 



and £4 is 



(1 -I- uuy 

_ ^ {d^udij^uf - {d''ud^u){d''ud^u) 
^~ {l + uuY ' ^ ' 

Further, m is a complex field which parametrizes the stereographic projec- 
tion of the target 5^, see Appendix A. The Faddeev-Niemi model is the S"^ 
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restriction of Skyrmc theory and so circumvents Derrick's theorem, because 
it consists of two terms such that their corresponding energies behave oppo- 
sitely under a scale transformation. The existence of (static) soliton solutions 
for the lowest Hopf indices has been confirmed by numerical calculations ([3] 

' [6]). 

There arc two more models which can be constructed from the two 
Lagrangian densities £2 and £4 separately by choosing appropriate (non- 
integer) powers of these Lagrangians such that the corresponding energies 
are scale invariant. For £4 the appropriate choice is 

£afz = -(A)^ (4) 

and for this model infinitely many analytic soliton solutions were found by 
Aratyn, Ferreira and Zimerman (=AFZ) by using an ansatz with toroidal 
coordinates ([7], [8]). We shall, therefore, refer to this model as the AFZ 
model in the sequel. The analysis of the AFZ model was carried further in 
([9]), where, among other results, all the space-time and (geometric) target 
space symmetries of the AFZ model were determined, and, further, the use of 
the ansatz with toroidal coordinates was related to the conformal symmetry 
of the model (more precisely, of the static equations of motion). It turns 
out that the AFZ model has infinitely many target space symmetries and, 
moreover, realizes the notion of classical integrability in a rather strict sense, 
because the static equations of motion (e.o.m.) resulting from the ansatz 
with toroidal coordinates may be solved by simple integration. Let us also 
mention here that the model with Lagrangian £4 has the same (infinitely 
many) target space symmetries and, in addition, a scale invariant action, 
leading thereby to analytic time-dependent solutions [10]. 
The other model is 

>CNi=(>C2)i (5) 

This model has first been proposed by Nicole ([11]), and it was shown in 
the same paper that the simplest Hopf map with Hopf index 1 is a soliton 
solution for this model. We, therefore, call this model the Nicole (=Ni) 
model in this paper. This model shares the conformal symmetry with the 
AFZ model, therefore, again the ansatz with toroidal coordinates may be 
used to simplify the static e.o.m. (to reduce them to an ordinary differential 
equation). However, the Nicole model only has the obvious symmetries - 
the conformal base space symmetries (in the static case) and the modular 
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target space symmetries, see [12]. Consequently, the e.o.m. are no longer 
integrable, and the solutions are no longer available in closed, analytic form 
(except for the simplest case with Hopf index one^). It is the main purpose 
of this paper to analyse these soliton solutions with higher Hopf index by 
developing the analytical treatment as far as possible and by performing 
numerical calculations where it is necessary. 

The AFZ and Nicole models currently do not have direct physical applica- 
tions, but they can serve as useful test labors for the understanding of generic 
features of models with knot solitons. This is all the more true as analj^ic 
results - already for solitons, but even more so for time-dependent problems 
like soliton scattering - are notoriously difficult to obtain for realistic mod- 
els like the Faddeev-Niemi model. Both the AFZ and the Nicole model are 
much easier to treat, and for time dependent problems the Nicole model even 
has one advantage over the AFZ model. The infinitely many symmetries of 
the former - which provide its integrability on the one hand - imply that the 
moduli space for multi-soliton configurations is infinite dimensional. This 
most likely makes an adiabatic treatment of soliton scattering problematic. 
On the other hand, that problem is absent for the Nicole model which may, 
therefore, probably serve as a natural test labor for an adiabatic study of 
the scattering of knot-like solitons. A further issue which may be useful for 
the study of time dependent problems is provided by the observation that all 
field configurations within the toroidal ansatz belong to an integrable sub- 
sector of the Nicole model (i.e., to a sector with infinitely many conservad 
currents), see Section 3. 

The paper is organized as follows. In Section 2 we provide the static e.o.m. 
of the Nicole model as well as the non-linear ODE resulting from the ansatz 
with toroidal coordinates. Then we analyse the latter equation in detail, 
establishing some properties analytically. In a next step we calculate the 
energies of the corresponding solitons numerically with high precision. We 
provide all energies for solitons with sufficiently small Hopf index. Further, 
we calculate some energies for larger Hopf index in order to establish some 
generic behaviour. We prove that this generic behaviour is in accordance 
with exact energy estimates for the ansatz in toroidal coordinates. 

In Section 3 we show that for all the solutions of Section 2 there exist 
infinitely many conserved currents, that is to say, all these field configurations 

^Exact solitons with higher Hopf index have been found in modified Nicole models [13] 
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belong to an intcgrablc subscctor of the Nicole model. 

In Section 4.1 we prove an upper bound for the energies of field configura- 

3 

tions with given Hopf index, of the type Eq < and provide an explicit 

value for the constant C2. In Section 4.2 we derive an analogous bound for 
the AFZ and Faddeev-Niemi model, which is easily achieved with the help 
of the results of Section 4.1. The discussion of the possibility to find a corre- 
sponding lower bound Eg > CiQ^ for the Nicole model, analogously to the 
Vakulenko-Kapitansky bound for the Faddeev-Niemi model, is presented in 
Appendix C. We show where a Vakulenko-Kapitansky type proof needs to 
be refined for the Nicole (and also for the AFZ) model and leave this lower 
bound conjecture for the moment. 

Section 5 contains our conclusions, where we comment on the relevance of 
our results for issues like stability, and on the relation to the corresponding 
results for the other two models (Faddeev-Niemi and AFZ). In Appendix A 
we collect useful results and facts about Hopf maps which we use throughout 
the paper. Appendix B contains some proofs which we need in the main 
sections. In Appendix C we discuss the lower energy bound. 

2 Solitons in the Nicole model 

2.1 Equation of motion 

The static e.o.m. for the Nicole model (5) is 

^{ujkUjUk + UjkUjUk) + UjUjUkk - {ukUkU + SukUku) = (6) 

where Uk = d^u etc., and the Einstein summation convention is understood. 
We now introduce toroidal coordinates (?], ^, ip) via 

X — g~^sinh?7C0S</7 , y — smh.r]sunp 

z — g~^sin^ ; q — cosh, rj — cos ^. (7) 

Further, we choose for u the ansatz 

u = f{ri)e''^'^+'''^ , m.n^TL (8) 

which is compatible with the e.o.m. for the Nicole model as well as for the 
AFZ model. The reason for this ansatz may, in fact, be understood from 
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the symmetries of the modeL More precisely, the presence of the conformal 
symmetry on base space imphes that the ansatz (8) is an "educated guess" 
for a solution to Eq. (6) in the sense of the Lie theory of symmetry. That is 
to say, if we choose a rotation about the z axis and a certain combination of 
proper conformal transformation along the z axis and translation along the z 
maximal set of two commuting base space transformations, then the 
corresponding infinitesimal symmetry generators (vectors v*) are precisely 
given by the tangent vectors along (/? and ^, = 9,^, and = 9^. The ansatz 
(8) is invariant under a combination of these base space transformations and 
phase transformations of the target space variable m, i.e., under the action 
of the vector fields = 9,^ — imudu and = (9^ — inudu, which provides 
precisely the educated guess according to Lie. A concise discussion of these 
points can be found in Ref . [9] , where the symmetries of the AFZ model are 
discussed in great detail. 

We need the gradient in toroidal coordinates 

V = {Vv)d, + {VOd^ + {Vip)d^ = q{e^d, + e^^^ + -^e^d^) (9) 

where (e^, e^, e<^) form an orthonormal frame in IR^. Further we need the 
relations 

1 — cosh 77 cos ^ 

V • e„ = -sinhTyH ^, V-e£ = -2sin^, V • = 0. (10) 

smhr/ 

Inserting now the ansatz (8) into the static e.o.m. (6) we find, after a straight- 
forward calculation, that the ansatz is indeed compatible with the static 
e.o.m. and that / has to obey the non-linear ODE 



2 

m 
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\ smh 7] J ^ ^ smh r] \ smh r\ ) 

where /' = 9^/ etc. It may be checked without difficulty that the simplest 
Hopf map 

u = sinh e*^"'""'' =^ / = sinhry , m = n = 1 (12) 
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indeed solves the above equation. 

For technical reasons it is preferable to introduce the new variable 

t = sinh?7 (13) 

in terms of which the above e.o.m. can be expressed as a pure polynomial in 
the independent and dependent variables t, f, ft, fu, 

F{t,f,ft,fu) = 

t'{l + t'){l + n {2t\l + t')f^ + {nh' + m')f) fu - U\l + effft+ 

+t\i + + + f^)f^ - 2t\i + e){nH'' + m'')f^f^+ 

+e p + n\l + 2^)) (1 + f')ff\ {nH^ + m^ff\\ ~ f') = (14) 

where ft = dtf etc. Further, we need the energy functional E[f] which results 
for the ansatz (8) after the integration w.r.t. the variables ^, (p. With 

dV = d^r — q'^sinhrjdrjd^dip (15) 

we find 



poo 

E[f] = 327r2 / dri sinhr] 
Jo 



For n = m = 1 for the solution f{t) = t the energy may be calculated 
analytically, 

If we calculate the energy for the field configuration f{t) — ct instead, we 
find 

E = 327i^V2-{c + -) (18) 
2 c 

which certainly has a minimum at c = 1 but also shows clearly that the one- 
soliton sector is indeed separated from the trivial sector (with Hopf index 
zero) by an infinite energy barrier. 



6 



Next we have to fix the boundary conditions which / has to obey. The 
u of Eq. (8) has to take values in the whole target space (Do in order to be a 
Hopf map. This implies that / has to take values in the whole positive real 
numbers, including zero, i.e., in the whole IR+o- It follows that the possible 
boundary conditions on / are restricted to one of the two following options. 
Either /(O) = 0, /(oo) = cxo or /(O) = cxo, f{oo) = 0. The reason for this 
is that the values / = and / = oo correspond to the north and south pole 
of the target space S'^, respectively. Consequently, their pre- images in the 
base space IRq must be one-dimensional objects (closed lines). But the only 
two values f] = rjo = const (or t = to = const) which correspond to lines 
rather than surfaces (tori) are the values r] = {t = 0) and rj = oo {t = oo). 
Therefore we assume 

/(^ = 0) = /(i = 0) = 0, /(r^ = oo) = /(i = oo) = oo (19) 

which is general, because the other option is related to this one by a symmetry 
transformation. 

2.2 Asymptotic behaviour and numerical evaluation 

In a next step we want to determine the asymptotic behaviour of the function 
f{t) for small and large values of t from the differential equation (14). For 
small t we assume that f{t) ~ r + o(i^) where (3 > a > 0. From Eq. (14) 
we get 

F{t,f,ft,ftt) ~ {2a''-a^+a^m^-am^-m'')t^"+o{t^''+^)+o{t^>^) = (20) 

The condition that the leading (i.e., smallest) order is absent therefore 
leads to 

2q;^ -a^ + a^rr? - aw? - = (21) 

which has four solutions for a of which only one is acceptable (real and 
positive), namely 

= ^ (l + V8m2 + l) (22) 

In the same fashion we may determine the subleading (higher order) contri- 
butions to / for small t in an iterative manner. We find 

/ = t"'"P(°)(t,t"'") for t«\ (23) 
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where P(°)(t, is a polynomial of its arguments, which is determined up 
to a multiplicative constant, 



p(o)(i,i"-) = ^ 4';^i¥"-. (24) 

fe,/=0 

Here it is possible to determine the higher coefficients c'^i in terms of m, n 
and Cqo^. Unfortunately, it is not possible to determine Coo\ i.e., an overall 
constant remains undetermined by this asymptotic analysis. 

We repeat the same asymptotic analysis for large values of t. We assume 
that f{t) ~ + o(i^) where now a > /? > and find 

F{t, f, ft, ftt) ~ (2a^ -a^ + aW - an^ - + o(i'"+') + o{t''^+^) = 

(25) 

and therefore the same condition for a like for small t with the only replace- 
ment m ^ n. We may again determine / asymptotically up to an overall 
constant, 

/ = i""P(°°)(r\r°^") for t»l (26) 

where 

oo 

p(°°)(r\r"") = cift-H-^'^- (27) 

k,l=0 

and again the higher coefficients c^^^ are determined in terms of m, n and 
c^Q. Again, Cq^^ remains undetermined. 

Next we turn to the numerical evaluation of the soliton energies for general 
m, n. For this purpose it is preferable to minimize the energy functional 
(16). The problem is that standard cvohition procedures for the differential 
equation (14) cannot be used because of the singular nature of this equation 
ait — 0. We proceed as follows with the minimisation. We first factorize the 
leading behaviour, e.g., we choose 

f^t)^t'^-{l + t')'^g{t) (28) 

where g obeys 

g>0 V t, g{0)=go, g{oo) = g^ (29) 
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Then we make a finite parameter ansatz for g of tlie form 



where 5^0; 9oo-i and hj are the parameters with respect to which we minimize 
the energy functional (16). For each value of m, n we increase the number 
of parameters (i.e., the integer number A;), until we reach stability (i.e., until 
the energies do not change further with the increase of the number of pa- 
rameters). We do not choose a full numerical minimization of the energy 
functional for the following reason. It turns out that the leading behaviour 
(i.e., g{t) = const), already provides rather good results for the correspond- 
ing energies. This implies that the energy functional, when viewed as a 
functional of g, is rather shallow. As a consequence, a full numerical min- 
imisation which implements the derivatives numerically has the tendency to 
produce false minima unless the numerical minimisation grid is chosen very 
finely (and consumes a lot of computer time). On the other hand, within our 
finite parameter ansatz we are able to perform the derivatives analytically, 
and therefore this problem is absent. 

We show our numerical results in Table 1 for all (m, n) = (1,1)... (5, 5) as 
well as for some higher values. Further we plot the energies versus the Hopf 
index Q = mn in Fig. 1. From our results it is obvious that the energies grow 
linearly in Q to a very good accuracy. More precisely, the energies for m = n 
lie on a straight line in Fig. 1 almost exactly. The energies for m 7^ n lie 
slightly above, whereas no energy lies below this hue. It is in fact possible to 
prove the linear growth of energies with the Hopf charge for solutions of the 
ansatz (8). This we show for the simplest case m = n in the next subsection. 

On the other hand, for the full model the sublinear upper bound E < 

3 

6*2(54 holds, which we prove in Section 4. It follows that the solutions we 
found for the symmetric ansatz (8) cannothe true solitons, i.e., global minima 

of topological sectors with fixed Q for sufficiently large values of Q. This 
is akin to the situation in the Faddccv-Niemi model and differs from the 
situation in the AFZ model (see our discussion in Section 5). 

2.3 Upper and lower bounds for the torus ansatz 

In this subsection we prove the bounds ciQ < E[fm,n] and, for the special 
case m = n, E[fm,m] < C2m^ (with explicit values for the constants ci and 
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n 


m 


^* 


n 


m 


E* 


1 


1 


1.41362193 


4 


1 


6.52450266 


1 


2 


2.62615091 


4 


2 


9.03384901 


1 


3 


4.13338723 


4 


3 


12.7436624 


1 


4 


5.37349409 


4 


4 


16.789796 


1 


5 


7.54150789 


4 


5 


20.886948 


2 


1 


2.64753548 


5 


1 


9.48680075 


2 


2 


4.52249973 


5 


2 


11.8901415 


2 


3 


6.62925943 


5 


3 


15.9786399 


2 


4 


8.85370929 


5 


4 


20.8937915 


2 


5 


10.9513184 


5 


5 


25.975942 


3 


1 


4.28232266 


6 


6 


37.2017869 


3 


2 


6.64941801 


8 


8 


65.771712 


3 


3 


9.64059435 


9 


9 


83.11675 


3 


4 


12.7307136 


10 


10 


102.501125 


3 


5 


15.8672836 


20 


20 


408.581244 








5 


10 


51.8990747 








6 


15 


96.0585933 








11 




103.670713 



Table 1: The rescaled energies E* = (327r^) ^E of solutions for selected 
values of m and n. 
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C2) in order to further support our numerical results. Here E[f] is the energy 
functional (16) and fm,n is the corresponding solution of the equation of 
motion. We first prove the lower bound. We use the inequality 

{a + b)'^>^ah for a,b>0 (31) 

which we prove in Appendix B (here is an optimal value). We choose 

a = and 6 = (n^ + vr? It^)p j (1 + t^) and introduce 

vn 

Q — mn , jjL = — (32) 

n 



to get 

E[f] = 327r2 / dtt{l + 1 



{f! + Q{i + §)li 



10 




> 327r^g3V3 / dtft—^^ = E[f] (33) 



where we used {t/ii + /i/t) > 2. Up to now we have in fact used inequalities 
for the energy densities, i.e., these inequalities hold for given, fixed functions 
/. However, if we evaluate the last expression for its minimizer f^^^ then 
the inequality holds for arbitrary / and, consequently, for the solutions fm,n- 
The minimization of the last expression is simplified by the observation that 
it is a total derivative. Therefore it is minimized for all / which obey the 
required boundary conditions /(O) = and /(oo) = 00. Explicitly we get 



/=oo 

E[f] = 327r2Q3V3 



/ / 1 ' 

H — — + - arctan / 



4(i + /2)2 8(1 + /2) 8 



/=o 



= 327r2g3\/3— = 1.02 X 327r2g. (34) 

Comparing with Table 1, we see that for larger values of m, n our lower 
bound gets very close to the numerical values, especially for m = n (e.g., 
better than one per cent for m = n = 10 or m = n = 20). This demonstrates 
the very good accuracy of our numerical results. 
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For the derivation of an upper bound we restrict to the simpler case 
m = n. We simply insert a trial function into the energy functional 

E[f]. Concretely we choose /^^^ = and get 

/•oo ™3^3m— 3o^ 

E[n^32n^l dtt{l^e)^^^-^ (35) 

and, with the substitution u — f^, 

Eit"^] = 327rW2i T du (u^ + m"^) , 

Jo ^ (1 + ) 

= 327rW2i-fl-^')^^-. (36) 

The function {1 — 1/ w?) / (cos 2^) is a monotonously decreasing function in m 
and may, therefore, be estimated by its value at m = 2 for the range m > 2. 
With 

(1 - ^ = (37) 

we therefore get 



4/ cos f 



< 327r2m2— = 1.18 x 327r^m^ for m > 2 (38) 

8 

which, for large m, is about 15 percent above the numerical value. By using 
Eq. (36) instead of the order rn? estimate (38), we get slightly better results 
which, for large m, are less than 10 percent above the numerical values. 

In short, we established the bounds CiW? < E[fm,m] ^ C2m^ with Ci = 
1.02 X 327r^ and C2 = 1.18 x 327r^, as announced. 



3 Infinitely many conserved currents 

In this section we want to demonstrate that the ansatz (8) belongs to a 
subsector of the Nicole model with infinitely many conserved currents for 
arbitrary profile function f{r)). The integrability condition defining this sub- 
sector exists for a wide class of models, and discussing it for the whole class 
does not complicate matters (for a more detailed discussion, and especially 
for the geometric meaning of the integrability condition, we refer to [14]; the 
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concept of higher dimensional integrabihty, which is at the basis of these dis- 
cussions, was introduced in [15]). Therefore, we start with the general class 
of Lagrangian densities 

£{u, u, M^, M^) = F{a, b, c) (39) 

where 

a — uu, b — Ufju'^ , c = {u^jU^Y — u^^ul (40) 

and F is an arbitrary real function of its arguments. Obviously, the Nicole 
model belongs to this class, = [4(1 + a)~^6]2. 
Further, we define the currents 

K^" = h{a)U^' (41) 

where h is an arbitrary given real function of its argument, and 11'* and Tl^ 
are the conjugate four- momenta of u and u, i.e., 

= £„M = u^Fb + 2(m\a«m - ulu^)Fc. (42) 

Finally, we define the infinitely many Noether currents 

= iiGuK^ - G^K^) (43) 

where G is an arbitrary real function of u and it, and G„ = d^G. 

Now we want to study the conditions which one has to impose in order 
to make the divergence df^J^ vanish. A simple calculation reveals 

d^J^ = i/i(^[(^txG'„ + G'„>J-(^iiG'^ + G'^^)iZj]F6 

+ (uG^ - uGu)[j{bF, + 2cFe) + Fj j (44) 

where the prime denotes the derivative with respect to a. 

The second term at the r.h.s. of Eq. (44) certainly vanishes if 

uGu -uGu^O (45) 

with the general solution 

G{u,u) =G{uu) = G{a). (46) 
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Without taking condition (46) into account, the first term on the r.h.s. of 
Eq. (44) vanishes for the integrability condition 

ul = (47) 

i.e., the complex eikonal equation (it also vanishes for some more complicated 
integrabihty conditions which we do not discuss here, see Refs. [16], [17]). 
However, by using condition (46) we may re-express the first term like 

{h'C + hG")Fb[u\l - u^ul] (48) 

and, therefore, we find, instead of the complex eikonal equation, the weaker 
integrability condition 

u\l-u^ul = Q. (49) 

The meaning of this condition becomes especially transparent when we re- 
express u in terms of its modulus and phase like 

w = exp(E + iA). (50) 

Then the complex eikonal equation is equivalent to the two real equations 

K - K (51) 

and 

E'^A^ = (52) 

whereas the weaker condition (49) becomes Eq. (52) alone or, for time- 
independent u, 

(VE) • (VA) = 0. (53) 

Remark: The condition G{u,u) = G{uu), Eq. (46), restricts the space 
of allowed G to a subspacc which is still infinite-dimensional, therefore the 
integrability condition (49) really defines a subsector with infinitely many 
conserved currents. 

Remark: The ansatz (8) obeys the integrability condition (53) for arbi- 
trary profile function f{rj), as is obvious from the orthonormality of the basis 
vectors {e^,e^,e^). Therefore, all our solutions of Section 2 really belong to 
the subsector of the Nicole model defined by condition (49) with infinitely 
many conserved currents. This is in contrast to the complex eikonal equation 
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(47), which leads to a non-hnear ODE for the profile function which is not 
compatible with the Nicole model field equation except for the simplest Hopf 
map m — n — 1, see [18]. 

Remark: Both condition (46) and the integrability condition (49) do not 
restrict the Lagrangian density, therefore the corresponding integrable sub- 
sectors with infinitely many conservation laws exist for all Lagrangians of 
the type (39). So they may be of interest for other models like, e.g., the 
Faddeev-Niemi model. 

4 Energy estimates for the full model 

In this section we study upper bounds on the energies of solitons with a given 
Hopf index, of the type Eq < C^Q"^ ■ In Subsection 4.1 we prove the upper 
bound Eq < C2Q^ by following the method used in [19] for the Faddeev- 
Niemi model. Further, we provide an explicit value for the constant C2. In 
Subsection 4.2, we briefly derive the analogous upper bounds for the AFZ 
and the Faddeev-Niemi models, again with explicit values for the constants 
C2. They easily follow from our results from Subsection 4.1. The issue of 
a lower bound CiQ^ < Eq is discussed in Appendix C. There we show 
where a naive Vakulenko-Kapitansky type proof does not work and needs 
some reflnement in our case, leaving open for the moment the problem of 
rigorously establishing such a lower bound. 

4.1 The estimate Eq < C^Q^'^ 

In this section we use the unit vector n to denote Hopf maps, because this 
has some technical advantages (components of the vector n can be easily 
estimated, whereas the complex function %i may take arbitrary values). In 
our proof we will closely follow the method used in [19] for the proof of an 
analogous result for the Faddeev-Niemi model. However, we will use explicit 
expressions for certain functions instead of using just their existence, which 
enables us to provide an explicit value for the constant C2 in our estimate. 

For a Hopf index which is a square, Q = P, a Hopf map is constructed 
in [19] which is a composition of a Hopf map n{f) with Hopf index 1 and 
a map N{n) : 5*^ — > S*^ with winding number w = I. The Hopf map n{f) 
is constructed such that it varies smoothly from n — (0, 0, 1) at r = to 
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n = (0, 0, —1) at \f\ = R and remains at this value for r > R. This imphes, 
by Eq. (114) of Appendix A, that it has indeed Hopf index 1. Further, the 
gradient |V^^fi| may be estimated hke jV'^^nl < c/ R ior r < R and by 
|V^^ri| — ioY r > R (here and below c denotes an unspecified constant). 
The map N{n) : S"^ ^ S'^ is constructed such that N varies smoothly over 
I nonintersecting geodesic discs covering the base S"^, where N = (0,0, 1) in 
the center of each disc and = (0, 0, —1) at the boundary of each disc and 
between the discs. Each disc contributes a winding number 1 for the map 
7V(n), therefore the total winding number is I. Further, the geodesic radius of 
each disc can be chosen (not larger than) Rg — d~2 in order to be able to put 
/ nonintersecting discs on one unit S*^. The gradient may be again estimated 
by the inverse (geodesic) radius, jV"^ A^(n)| < c/2, and with the help of the 
chain rule one finds for the gradient of A^ that \'V^oN{n{r))\ < cR~^l^ and 
for the energy density £^Ni(^) ~ |V^oiV|3 < cR-H'-^. The energy Ej^i, which is 
the integral over all space of Emir), may in the case at hand be calculated by 
integrating over the ball of radius R with the result that E^i < = cQ^ 
which is the announced estimate for Q = 1"^ . 

Before rederiving this result in more detail and with explicit choices for 
the maps described above in rather general terms, it may be useful to have 
a geometric picture of the above compostion of maps. In fact, the map 
N{n{r)) describes I nonintersecting, full tori confined to the ball B^, where, 
in addition, each full torus is linked with all other tori. In the center of 
each torus there is a closed line along which A^ takes the value A^ = (0, 0, 1). 
From the center to the surface of each full torus, A^ varies smoothly from 
TV = (0, 0, 1) to iV = (0, 0, -1). At the surface of each full torus (which has 
the topology of the two-torus T^) as well as in the space between the tori 
and outside the ball B^, N remains at the value A^ = (0, 0, —1). 

Now for the detailed derivation, where we still assume Q = P for the 
moment. For the map n{r) = (sin -(9 cos 0, sin -i? sin 0, cos -i?) with Hopf index 
1 we assume that it has cylindrical symmetry, as in Eq. (118) of Appendix 
A. For the profile function p(r) we want to find a choice which varies as 
smoothly as possible from to tt as r varies from to R, in order to be able 
to estimate it with a number as small as possible. The smoothest choice 
p = nr/R for r < R and p = n for r > R does not have a continuous first 
derivative at r = i?, which we need for the energy density, but there is a 
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simple generalization which does have one. We choose 



p,(r) = 7rF,{x) , x = ^ (54) 



where 



2 

F,{x) = F<(x) = - — -X for < X < 1 - e 

F,(x) = F>{x) = 1- -^^^^—^{x-lf for l-e<x<l 

F,{x) = 1 for x>l (55) 

and e is a sufficiently small, positive nonzero parameter. It may be checked 
easily that Pc{r) is continuous and has continuous first derivatives everywhere. 
Further, the following estimates hold, 

\F,{x)\ < k{e)x for 0<x<l 
\Fe{x)\ < k{e) for 1 < x < oo 
\F^{x)\ < k{e) for < x < oo (56) 



where 



For the map N{n) with winding number / we choose / nonintersecting 
geodesic discs on the base S'^ such that varies smoothly from = (0, 0, 1) 
in the center of each disc to N = (0, 0, —1) at the boundary of each disc and 
remains at this value in between the discs. The geodesic radius ■j?o of these 
geodesic discs has to be chosen sufficiently small so that it is possible to put 
/ nonintersecting discs with this geodesic radius on the base S^. We prove in 
Appendix B that the choice 

*o = j^ (58) 

— * 

is sufficient. Next we give an explicit expression for the map N for one 
geodesic disc which is chosen symmetric about the north pole. With 

N — (sin X cos cr, sin x sin a, cos x) (59) 

and 

n — (sin 1? cos 4>, sin sin 0, cos i?) (60) 
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we get 

xW=7rF,(^), ^ = ^ (61) 

and 

(7 = (62) 

where F^{-) is defined in (55) (with the additional restriction that here i? < 
Ay/l because < n). The above-defined N has winding number 1. and 
to reach winding number / we have to require that it is nonconstant over 
/ — 1 further geodesic discs. We do not need, however, explicit expressions 

— * 

for the contribution of these further discs to N, because their gradients on 
S'^, \W^'^N\ can be estimated by the same majorants as the above exphcit 
expression. The reason for this is that the above symmetrically chosen disc is 
related to an arbitrary disc with the same geodesic radius by a rotation of the 
base space S^, and the gradient | V*^ is invariant under such rotations. All 
that changes is that more regions of the base S"^ provide nonzero contributions 
to the gradient. 

Next we want to estimate the energy density S-^i{r) ~ |ViVp. Here we 
first estimate (Einstein summation convention is understood, k,a — 1 ... 3) 

I VTVp = iV.mf = xU'^k^^f + sin^ x(V,a)2 (53) 

where we used (59), and = d^x- The first term may be estimated by 

Xl<{l)'k{er (64) 
and, with the help of spherical polar coordinates, by 



Apf cos^ p sin^ 9 ^ sin^ p 4 cos^ 9 



1 - sin^ p sin^ 9 1 - sin^ p sin^ 9 

< 4p? + 4^<8vr^^ (65) 



where we used 



cos^ a 



1 — sin a sin /3 



2 < 1 , sin a < a (66) 



18 



as well as Eqs. (54) and (56). For the second term we find 



sin^x < sin^ Uk{e)pj < ('^V sin^i? (67) 



where we used 



TT 

sin ax < asinx for a > 1 , x E [0, — ] (68) 



which we prove in Appendix B. Further we find, using a = (j) and the expres- 
sion for (p in Eq. (118) 

sin^^?(Vfe0)^ = 4sinVsin2^(l -sinVsin^^*)- (69) 

/ 1 p^cos^^ + (l/r^)cos^psin^psin^^\ 
■ Vr^sin^^ ^ (l-sinVsin'^)2 ) ^ ' 

^ 4^(1- sin Vsin^^) (71) 
cos^ 

+4:sm^ psin^Opl- (72) 

1 - sm"' p sm^ ^ 

, . • 2 • 4/)SinV cosV 

+4sin psm — ^ 0-77 (73 

r2 l-sinVsin^^ ^ ^ 

Putting everything together, we therefore find 

Taking the power of | and integrating over the ball we estimate the 
energy by 

E^, = / d'r\WNf < [ dhhk^ { = "^^lil^si^kieff (76) 
J Jb^ \ R I 3 
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Finally, for the estimate it is sufficient if the estimating field configuration is 
a distribution. Therefore we may now perform the limit lime^o k{e) = 1 and 
get, for Q ^f, 

Em < y (87r)35i/i (77) 

In the cases when the Hopf index Q is not a square, we shall find a 
slightly weaker bound (i.e., a slightly larger value for the constant C2). Again 
following [19], we write Q as 

Q = f + m, me [1,21] (78) 

and estimate the energy by the following Hopf map. It maps the ball 
with radius R about the origin to the two-sphere with a contribution of 
to the Hopf index in exactly the way constructed above for the case Q = P. 
Further, it maps m balls Bl, i = 1. . .m with unit radius to the S"^ such 
that each ball contributes one unit to the Hopf index. Here the m unit balls 
are chosen such that they intersect neither each other nor the ball B^. The 
map from B^ contributes exactly the above-calculated expression (77) to 
the energy, whereas the m maps from the Bj contribute the energy of a unit 
Hopf map each. One might think to take just the energy (77) with / = 1 
for each energy contribution, but we can do better. The reason is that we 
were forced to choose ■j?o rather small in Eq. (58) to ensure that we can put 
/ geodesic discs on one S^. But here we are only interested in a map with 
winding number one, therefore we can choose the "geodesic disc" equal to 
the full S'^, i.e., t^o = ti"- This leads to the following energy estimate 

E^, < —(87r)^5kl + — (27r)^5im (79) 
3 3 

which we want to estimate in terms of Q again, i.e., 

C2/^ + Cam < C4(/^ + m)^. (80) 

Here the l.h.s. grows linearly in m, whereas the r.h.s. grows sublinearly, 
therefore the inequality certainly holds for all m € [1,2/] if it holds for the 
maximum value m = 21, i.e., 

C2l^ + 2C3I < 04(1'^ + 2l)i . (81) 
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With the help of the inequahty 

1 2 

C2 + 20^1-^ < C2 + 2C3 < (C2 + 2C3)(1 + y)^ (82) 

we see that the above inequahty certainly holds for C4 — C2 + 2c3. Inserting 
the concrete numbers from (79) we therefore get 

^Ni<y(27r)3(4=^ + 2)5igt (83) 
as our final estimate for general Q. 

4.2 Upper bounds for the AFZ and Faddeev-Niemi 
models 

We may use our result to obtain upper bounds for the energies of the AFZ 
and Faddeev-Niemi models, as well. For the bound on the AFZ model we 
use the simple inequality for the energy densities 

£4 < (84) 

which is obvious from the ('(9,0) parametrization, see Eqs. (119) and (120). 

3 

It follows immediately that Eafz < (|) * -^Ni and, therefore, 

3 

Eafz < (^)' ^^^Trfii' + 2)5igt. (85) 

For the Faddeev-Niemi model with energy density £^fn — £2 + we 
consider the case Q = of a square Hopf index first. We use Eq. (75) for 
the density S2 and the inequality (84) and find 

Efn = / d\i£2 + \£4) 




and, with the choice 

R = RoVl (87) 



we get 



^^FN < 20|^(47r)2zi fi?o + 1° W;^) ■ (88) 



For a neat estimate we now minimize the expression in brackets w.r.t. Rq 
which leads to 

i?o = 47rv^ (89) 

and to 

^FN < 40^(47r)^yiOAZi. (90) 

Here we have separated a factor (47r)^ because it can be replaced by tt^ for 
winding number one, exactly as in the estimate for the Nicole model, which 
we need for general Q. Again we write Q — P -\-m and estimate 

47r , , , n , — -, 3 , An 



£^FN < 40— (47r)^\/lOAZ5 + 40— Tr^vToAm (91) 



and finally by 

An , 



Efn < 40— 7r-^V10A(4^ + 2)Q^ (92) 



in complete analogy to the Nicole model. 



5 Conclusions 

Combining numerical and analytical techniques, we found the energies for 
static solutions of the Nicole model within the symmetric ansatz (8) for a 
wide (in principle arbitrary) range of the integers m and n. It turns out 
that the energies of these symmetric solutions grow linearly with the Hopf 
index Q = mn within a very good accuracy. Together with the sublincar 
growth inequality of Section A, E < 6*2(5^ this implies that, for sufficiently 
large values of Q, the symmetric solutions we found cannot be global minima 
within their respective topological sectors with fixed Hopf index Q. It is of 
some interest to compare this result with the corresponding results for the 
AFZ and Faddeev-Niemi models. 
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For the AFZ model all solutions for the symmetric ansatz (8) and their 
corresponding energies can be calculated exactly [8]. The energies arc 



and, for m = n (/X = 1) perfectly fit within the upper bound £^afz < C2Q^ 
which holds also for the AFZ model (see Section 4). It is, therefore, plausible 
to conjecture that these solutions might be true minima for m = n (and, 
maybe, even for m ^ n provided that m and n do not differ too much). 

For the Faddeev-Niemi model the ansatz (8) is not compatible with the 
equations of motion (due to the lack of conformal symmetry of the latter), 
but the cylindrically symmetric ansatz 



for complex / is compatible and leads to an equation in the two independent 
variables r] and ^. Numerical investigation of the full e.o.m. in three vari- 
ables shows that the solutions for the ansatz (94) are true minima only for 
the lowest values Q = 1,2 of the Hopf index [4]. In this respect, therefore, 
the Nicole model bears more similarity with the Faddeev-Niemi model than 
with the AFZ model. It is tempting to speculate at this instant that the dif- 
ferent behaviour of the AFZ model - i.e. the possibility that the symmetric 
solutions (93) are true minima for arbitrary Q - is related to the integrability 
properties of the latter. But at the moment this is, of course, only a specu- 
lation which deserves further investigation. 

Appendix A: Facts about Hopf maps 

A Hopf map is a map from the three-sphere to the two-sphere S'^. The 
third homotopy group of the two-sphere is nontrivial, 113(5'^) = 2, therefore 
Hopf maps are characterised by an integer topological index, the so-called 
Hopf index Q. The three-sphere is topologically equivalent to one-point com- 
pactificd three-dimensional Euclidean space IRq, therefore, for each Hopf map 
— s> S"^ there exists a corresponding map IRq S"^, which we shall call Hopf 




(93) 



3 



(94) 
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map, as well, and which has the same Hopf index (due to the metric inde- 
pendence of a topological index). Explicitly, a Hopf map IRq S"^ may be 
given by a three-component unit vector field 



n(r) : IRq S'^ , = 1 , lim = no = const (95) 

f— +00 



where the tip of the unit vector field spans the unit two-sphere, or via stere- 
ographic projection 

^ f ^- -f -\ ^ -\ ni+in2 , . 

n = ; (u + u, —liu — u),l — uu) : u = . (901 

1+ I |2 ^ l + ns ^ ^ 

by a complex field 

u{r) : ]Rq — > (Co , lim = uq = const (97) 

f— »oo 

where our conventions are such that the projection is from the south pole to 
the equatorial plane of the two-sphere. Further, (Dq is the one-point compact- 
ified complex plane. A third possibility to parametrize a Hopf map is by the 
spherical angles G [0, 27r], 'd e [0, tt], which are related to the unit vector n 
and complex field u via 

■& 

u — iajci—e^^^ n = (sini?cos0, sin'j^sin^, cos'j?). (98) 

The geometry behind the Hopf map may be visualised as follows. The 
pre-images under the inverse of the Hopf map u of points in the target space 
S'^ (or, equivalently, (C^) are closed curves in IRq (in general, knots), where 
each two curves corresponding to two different points in target space are 
linked exactly Q times. An analytic expression for the Hopf index Q is 



^^Te^^y^^^-^ ^^^^ 

— * 

where B is the Hopf curvature 

B = 2 VlO^Vw ^ _\^ahc^a^^h X = - sini^Vi? X V0 (100) 

and A is the gauge potential for the "magnetic field" i3, i3 = V x There 
is no local expression for A in terms of the Hopf map u alone (however, there 
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certainly exist non-local expressions for certain gauges, as is obvious from 
the analogy with electrodynamics, by choosing, e.g., the Couloumb gauge 
V • ^ = which implies ^ = V x Jd^r'A-\f- r')^{r')). 

— * 

Geometrically, B is the Hodge dual of the puUback of the area two-form 
on S"^ under the Hopf map i.e., 

B = Bkdr^ (101) 
^ = ^^jkdr^dr'' , B^ = Skim^im (102) 

where 

The second cohomology group of the three-sphere is trivial, H2{S^) — 0, 
therefore the closed two-form must be exact, JF = dA. It follows that the 
gauge potential A is globally well-defined for appropriately chosen gauges 
(e.g., in the Couloumb gauge). 

The simplest Hopf map with Hopf index 1 is 

u = i 2(^ + i?/) 

2z + i{r'^-l) ^ ' 

(the irrelevant factor i has been introduced to be in exact agreement with 
the simplest Hopf map of Eq. (12)). Its level curves (i.e., the pre-images of 
points u = const) are circles which lie on tori, and each two different circles 
are linked precisely once. 

There are different ways to construct Hopf maps explicitly. One method 
which we shall need is to compose a given Hopf map with a map ^ S^, 

N{n{r)) : IR^ ""-^ ^-^^ (105) 

— * 

If the Hopf map n{r) has Hopf index Q and the map N{H) has winding 
number w, then the composition map 7V(n(r)) has Hopf index Q' = w^Q. 

A method to construct all possible Hopf maps starts from maps — > 
(or, equivalently, maps IRg — ^ 5*^). The third homotopy group of is 
nontrivial, n3(5'^) = 2, therefore such maps are classified by a topological 
index, the winding number W. It is possible to construct Hopf maps IRq — > 
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S"^ from maps IRg such that the Hopf index equals the winding number. 

Exphcitly a map IRq may be given by a four component unit vector 

field 

(ei(f),e2(f),e3(f),e4(f)) : 1R^^5^ e«e« = 1 , a = l,...,4 (106) 

lim Ba — e'^ — const (107) 

f— »oo 

then the corresponding Hopf map is given in terms of u as 

^.= ^i±^ (108) 

or in terms of n as 

fi = (2eie3 + 26264, -26i64 + 26263, 63 + 64 - 61 - 62). (109) 

Further, it is now possible to give an explicit, local expression for the gauge 

— * 

potential A, in terms of the 6^, as 
2 

^ = -[(ei - ^e2)V(6i + ^62) + (63 - i64)V(63 + ^64)]. (110) 

It may be checked without difficulty that indeed V x A = B where B is the 
Hopf curvature (100). 

As the group manifold of SU{2) is equivalent to the three-sphere S^, we 
may use maps IRq SU{2) instead. Indeed, with the group element 

U — 64 — ie- a = exp[—ipk ■ a] , lim U — Uq — const (HI) 

f—*oo 

64 = cosp, e — smpk, k —1 (112) 

(where e — (61, 62, 63) and B are the Pauli matrices) the Hopf map (109) may 
be expressed like 

ltr(73C/^aC/. (113) 

This latter representation may be used to produce ansaetze for Hopf maps 
with cylindrical symmetry from maps IRq — > SU{2) with rotational symme- 
try. Indeed, for 

k^r = f/r, p = p{r), p(0) = , lim p = ttW" , W E 'Z. (114) 
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(the conditions on p(r) are to ensure a well-defined U on all IRo, and they pro- 
vide, at the same time, a winding number equal to W for the corresponding 
map ml or ml SU{2)), the resulting group element U {p{r) , f) is ro- 

tationally invariant in the sense that any rotation in IRq can be compensated 
by a SU (2) transformation, 

U (p(r). Of) = cos p(r) — i sm{p{r)a ■ (Of) = V'^U {p{r), f)V (115) 

where O is a rotation matrix acting on f and V is the SU (2) matrix corre- 
sponding to the rotation O. Within this ansatz, and using spherical polar 
coordinates f = (sin 9 cos ip, sin 9 sin ip, cos 9), we get for the complex field u 

sin p sin ^e*"^ 

U^- ^ : 116 

sm p cos 9 + z cos p 

which has cylindrical symmetry in the sense that any rotation (p ^ ip + ipo 
can be compensated by a phase transformation u e*'^°M. Further, -u is a 
Hopf map with Hopf index equal to the winding number W. For technical 
reasons we prefer to use a u which is multiplied by —i in our main calculation, 
i.e., 

^ ^ + ^ sin p sin 9 gj(p+iarctan[tanpcos6»l (117) 

With u = tan(i?/2)e''^ this leads to 



sin = 2 sin psm9\Jl — sin^ p sin^ ^ , = (/? -|- arctan[tanpcos^]. (118) 
Finally, we display the energy densities of the Nicole and AFZ models in the 

3 

three different parametrizations, £m = £2 with 

£2 = \Vnf = I Vfenl^ = = (V^)^ + sin^ ^(V0)^ (119) 

3 

and Safz — £4 with 
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' (1 + uu)^ 



= sin' -d ({V'df{V(f)f - (Vi? • V<f)f) . (120) 



Appendix B 
Bl: Proof that 



3 1 3 

{a + b)^>cah for a, 6 > , c < V3-. (121) 



By squaring the above expression and subtracting the r.h.s. we get 

r' 

+ 30^6 + 3(1 - ci')a62 + ^3 > , = - (122) 

3 

or, by introducing the variable 3; = |, 

h{x) = + 3a;2 + 3(1 - d^)x + 1 > for x > 0. (123) 

This problem we solve by first calculating the minimum Xmin of h{x) for 
X > 0, and by then requiring that /i(xinin) = 0. The minimum is at 

/i'(x) = 3(x' + 2x + (l-d')) = =^ Xini„ = rf-1 (124) 

and the condition /i(xmin) = finally leads to 

/^(a^min) = (i^(3 -2ci) = A d>\ d^^ (125) 

which implies (121). QED 

B2: Proof that one may cover the unit sphere with I noninter- 
secting discs with geodesic radius = 

First we study the case where I is a square, I — k^. We divide the northern 
hemisphere of the two-sphere into segments by circles which are concentric 
about the north pole and are given by 

77777- 

m^l...k. (126) 
28 



This divides the northern hemisphere into one disc with geodesic radius ^ 
about the north pole (on which one geodesic disc with geodesic radius "f^o = 
certainly does fit) and into sphere segments such that the boundary circles 
have a geodesic distance of Therefore discs with geodesic diameter 2'j?o — 
^ exactly fit in between. Further, 2m + 1 geodesic discs fit on the segment 
between the circles at '(9 = ^ and at i!) = ^"^2k^^ ■ This is because the central 
circle of the segment has a circumference wich is bigger than the sum of m + 1 
diameters of the geodesic discs, i.e., 

2.sm ((„ + 1)^) > 2.^ [(m + i)iL) = (2„ + 1)^ > (2,„ + l)iL 

(127) 

where we used 

2 n 
smx>-x for xefO, -]. (128) 

TT 2 

Finally, summing over all discs on all segments (including the first disc at 
the segment which is itself a disc) we find 

k-l k 

53(2m + 1) - ^(2m- 1) = (129) 

m=0 m=l 

therefore we may indeed distribute k"^ nonintersecting geodesic discs with 
geodesic radius i?o = on the northern hemisphere. 

It remains to study the case where I is not a square. Then we write 
/ = k!^ + j. We distribute the k"^ discs on the northern hemisphere as just 
described. Further, we distribute the remaining j discs on the southern hemi- 
sphere. This is obviously possible as long as /c^ > j. So we cover all cases 
except the case / = 3. This latter case may be proven easily with the help of 
elementary geometry. QED 



B3: Proof that 



TT 

sin ax < a sin a; for a > 1 , x e [0, — ] (130) 

2 



We introduce a=l-|-e, e>0 and write 

sin ((1 + e)x) = sin a; cos ea; + cos a; sin ex 
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< sin X + ex cos x < sin x + e sin x = a sin x (131) 
where we used cosx < 1, sinx < x for a; > 0, and 

71" 

xcosx < sinx ^ tan a; > x for x e [0, — ] (132) 
which is obviously true. QED 



Appendix C: 

3 

On the possibility of a lower bound C\Q~^ < Eq 

It may be useful to briefly comment on the problem of deriving a lower 
bound for the energy of the type CiQi < Eq, analogously to the case of 
the Faddeev-Niemi model, where this bound can be proven (Vakulenko- 
Kapitansky bound, see e.g. [20], [21], [22]). In the case of the Faddeev-Niemi 
model, the first step in the proof is Hoelder's inequality (we use the conven- 
tions of [23], where all inequalities of this appendix can be found; further, we 
assume that A and B belong to the appropriate Sobolev spaces such that all 
integrals below exist) 

IGTr'g = J d'^rA- < (/ d'^r\A\''^ ' (^j ' (133) 

- + - = 1. (134) 

p q 

For the Faddeev-Niemi model one has to choose p — 6, q — ^ and may 
estimate the first term with the help of the Gagliardo-Nirenberg-Sobolev 
(GNS) inequahty, 

1 1 
(1 d^rlA^y <c(^J d'rMAW^y (135) 

(here and below c is an unspecified constant) and, when the Couloumb gauge 

— * 

condition V • ^ = is imposed, the integrand can be reexpressed like 

|V|^ir < c\^djAk\^ = c(|V X + V • A) (136) 

jk 
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where 

Ai = {A2d2 + ^3^3)^1 - ^1(^2^2 + 93^3) (137) 

and A2, A3 follow by permutation. Therefore, the integrand is the sum of 
and a total divergence which does not contribute to the integral (for A which 
decay sufficiently fast in the limit |f| — > 00; this is automatically satisfied by 

— * 

the class of A we consider). The proof continues with the interpolation 
inequality 



(/■ 



for the second term. This leads to 

167r'g < c( I dMB^^ 



dh\B\ 



£r\B\ 



(138) 



(139) 



— * „ — * 

and, together with \B\ — £4 and \B\ < ;^^2, to 



< cA~3 



A / d^rS. + / d^rSo 



d\S, 



(140) 



where we used Cauchy's inequality in the last step. The Vakulenko-Kapitansky 
bound follows immediately by taking the above inequality to the power |. 

Now let us point out where the analogous proof fails for the AFZ and 
Nicole models. One again starts with Hoelder's inequahty (133), but now 
one would have to choose p = 3 and q — With the GNS inequality one 
can again estimate the first term 



d\-\A\'' 



< c 



d'r\V\A\r^ 



2 

3 \ 3 



(141) 



but now the integrand is 



|V|^||i < c|^a,-A|^ = c(|V x^|2 + V-A) 

jk 



(142) 



and therefore no longer a total derivative, and we cannot express the in- 
tegrand in terms of the Hopf curvature alone. We have not been able to 
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overcome this difficulty. For our main results, however, the upper bound is 
more important, which shows that for sufficiently high Hopf index our solu- 
tions of Section 2 cannot be true minimizers of the energy. 
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Figure 1: The rescaled energies E* = (327r^)~^i? (vertical axis) are plotted 
versus the Hopf index Q = mn. The straight line connects energies for 
171 = 12. The energies for m ^ n lie slightly above this line. 
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